Thermally isolated Luttinger liquids with noisy Hamiltonians 
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We study the dynamics of a quantum-coherent thermally isolated Luttinger liquid with noisy Lut- 
tinger parameter. We calculate the out-of-equilibrium energy fluctuations as well as the correlation 
functions of this noise-driven system, and find universal features. By identifying two types of energy 
moments, which can help tease apart the effects of classical and quantum sources of fluctuations, 
we argue that the classical source is dominant in the thermodynamic limit. 
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Recent experimental developments with ultracold 
atoms have motivated numerous studies of the noncqui- 
librium dynamics of thermally isolated many-body quan- 
tum systems. Most of these studies focus on deter- 
ministic quantum evolution, generated, e.g., by a sud- 
den quench or gradual ramping of the Hamiltonian (see 
Ref. |l| and the references therein). However, stochastic 
driving of thermally isolated systems with noisy Hamil- 
tonians is much less studied, and, in particular, the 
role of quantum coherence remains largely unexplored. 
As shown in Fig. [TJ fluctuations in such systems stem 
from two distinct sources: i) the (classical) stochastic na- 
ture of the driving (different realizations of noise result 
in different wave functions with different energy expecta- 
tion values) and ii) inherent quantum fluctuations (each 
wave function can be a coherent superposition of states 
with different energies). 

The general setup of the problem is as follow. We con- 
sider a quantum system with Hamiltonian Hq = J2i9iOi> 
where <?iS represent the coupling constants and OjS are 
some local operators. The system is initially at zero tem- 
perature, i.e., in the ground state of Hq, where the initial 
energy distribution (see Fig. [I) is a delta function. For 
t > 0, the coupling constants fluctuate in time: gi{t) = 
9i+^9i(t)-, where, by assumption, \Sgi(t)\ <C \gi\. For each 
realization Sgi(t) of noise, the system is then described by 
a pure state, i.e., a coherent superposition, which evolves 
unitarily with Hamiltonian H(t) = Hq + J^i 5gi{t)Oi. 

Such noisy systems can be experimentally realized by 
cold atoms in optical potentials, with 5gi(t) generated by 
small fluctuations in the potential. For each realization of 
noise, the atomic system sees an effective time-dependent 
Hamiltonian H(t), which does not gencrically commute 
with Hq. Thus, the quantum evolution creates excita- 
tions with respect to the ground state of Hq, increasing 
the system's energy. Note that the driving with H{t) is 
assumed external, i.e., there is no feedback action from 
the quantum system on Sgi(t). Therefore, in isolation 
from a thermal environment (without a mechanism for 
dissipation), the system can absorb energy ad infinitum. 
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FIG. 1: (color online). The time evolution of (e) for two 
different realizations, I and II, of noise. Each (e) is the average 
of a quantum energy distribution. At t = 0, the system is in 
the ground state for all realizations of noise, a) Neither of 
the classical or quantum fluctuations dominate, b) Classical 
fluctuations dominate. 



The absorbed energy of such noise-driven thermally 
isolated systems is a time- dependent random variable e 
(as opposed to the steady states which emerge in noise- 
driven systems coupled to a heat bath [j|). What is the 
precise meaning of e for systems described by coherent 
superpositions of energy eigenstates? How can we char- 
acterize the average E(e) and the variance Var(e) of this 
random variable as a function of time? How are the cor- 
relation functions affected by this stochastic driving? 

Our objective in this paper is to address the above 
questions. To do so, we first identify two different types 
of quantum- and noise-averaged moments, which make 
the notion of a random absorbed energy precise, and, 
help tease apart the effects of the classical and quan- 
tum sources of fluctuations discussed above (see Fig. [T]) . 
We then perform explicit analytical calculations of these 
moments, as well as different noise-averaged correlation 
functions, for generic one-dimensional systems described 
by the Luttinger-liquid (LL) theory. 

Our main results are as follows. We find that, in the 
thermodynamic limit, the classical source dominates the 
energy fluctuations. In this limit, and in the regime of 
validity of the LL description, we find a universal rela- 
tionship (independent of the strength of noise and the 
Luttinger parameter) between the average and the vari- 
ance of energy: 



Var(e) = F(irt) [E(e)f 
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where T{x) is a dimensionlcss function of x = ntu/a 
(with velocity u and lattice spacing a set to unity), which 
decays as 1/x for large x. We also find that the noise- 
induced correction to equal-time correlation functions of 



scaling dimension A, C^°\x) 
universal form: 
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where W 2 is a constant with dimension of time [see 
Eq. ([8]) below] , which characterizes the strength of noise. 

Let us now discuss the two types of aforementioned en- 
ergy moments. As the evolution with different noise real- 
izations can lead to different pure-state wave functions, it 
is helpful to first introduce the notion of a wave-function 
probability distribution f(jp,t). We can simply think of 
/(?/>, t) as the joint probability distribution for all the in- 
dependent real numbers required to label \ip) modulo an 
overall U(l) phase. Similarly, J dip is shorthand for inte- 
gration over all such numbers. In the scenario considered 
here (the system initially in the ground state), f(ip,t) is 
a delta function at t = 0, which broadens for t > 0. Note 
that the wave functions i\) are, generically, coherent su- 
perpositions of all energy eigenstates: \tp) = J2 n c n\ n )i 
where Hq\h) = e rl |?i). 

The fluctuations of measured energies (with respect 
H ) are then characterized by the following moments: 
( e " i ) — J2n e ri'^'( e n), where the brackets (ovcrline) 
indicate a quantum (noise) average, and V(e n ) = 
j dij)f(ip,t)\ct\ 2 is the probability of measuring e„. This 
yields 



(e™> = J <h/>f(1>, *)(VWIV>> = tr [HS l p(t)] , (3) 

where p(t) = J di^f(t/j,t)\xp){ip\ is the density matrix at 
time t. The fluctuations described by Eq. ([3]) mix both 
classical and quantum contributions, which is also the 
case in experiments. 

However, we can alternatively consider the fluctuations 
of the expectation value of energy, (e) = (ip\Ho\^), over 
different realizations of noise, which are characterized by 



#/ty,t)(W>IW» r 



(4) 



As (e) is a quantum-averaged quantity, its fluctuations 
stem solely from the classical stochastic driving. Thus, 
the difference between these two types of moments can 
serve as a theoretical diagnostic for the relative impor- 
tance of the two sources of fluctuations (see Fig. [1} . Note 
that the moments above [Eq. (QJ] can not be written in 
terms of the density matrix for m > 1 @, 0] ■ Also notice 
that fluctuations of quantities other than energy can be 
similarly studied by expanding the wave function in the 
corresponding eigenbasis. 

Let us now turn to the the specific model studied in this 
paper, namely, a Luttinger Liquid, which is a universal 



low-energy description of interacting fcrmions and bosons 
in one dimension. Since a noise-driven dissipationlcss sys- 
tem can keep absorbing energy, this low-energy descrip- 
tion will eventually break down for most experimentally 
realistic scenarios. In this paper, we focus on dynam- 
ics over a finite time scale where the LL description re- 
mains valid. Notice that such time scales can be extended 
by decreasing the strength of noise. Luttinger liquids 
with negligible coupling to the environment, which sus- 
tain quantum coherence in their life time, have been ex- 
perimentally realized with cold atoms [8- 

In terms the Luttinger parameter K and velocity u, 
the LL Hamiltonian is given by 



h{k) = « W# n 9 n_ 9 + i q 2 

9>0 ^ 



(5) 



where $ g arc bosonic fields and Tl q their conjugate 
momenta. The Hamiltonian above can be written as 
J2 q >o + Hq^j , where Hf is the Hamiltonian of a 
single harmonic oscillator involving only the real (imagi- 
nary) part of $ 9 . It is convenient to shift the Hamiltonian 

by a constant, i.e., — > ^ — uq/2, so that the 
energies are measured with respect to the ground state. 
We consider a system initially in the ground state of Hq = 
H(K ), which evolves with H(K(t)) = H(K + SK(t)) 
for t > 0, where SK(t) <C K represents the noise. As the 
fluctuations of velocity u correspond to a trivial rescal- 
ing of the Hamiltonian, we set u to unity throughout this 
paper. 

Expanding the Hamiltonian ([5j in SK results in 
quadratic (in the bosonic fields) noise terms. Thus, inte- 
grating out the noise at the outset (as in Refs. [3I IR Il2j) 
leads to an interacting (quartic) effective action, which 
is difficult to treat exactly. We thus take an alternative 
approach described below. As the time-dependent Lut- 
tinger parameter is assumed spatially uniform, momen- 
tum is a good quantum number, and, the ground state 
of Hq is a direct product of Gaussian wave functions. 
Since each mode evolves with a quadratic Hamiltonian, 
the time-dependent wave function of the system, i.e., the 
solution of the time-dependent Schrodinger equation, re- 
tains the form 



=n 



2 9 [»«,(*)] 



CXp (-<? Z q (t) |$g| 



(6) 

where the parameter z q (t) now satisfies the Riccati equa- 
tion iz q (t) = -j^jj [(K(t) z q (t)) 2 — l], with initial condi- 
tion z q (0) = Kq 1 0, Hi- Note that each single-mode 
wave function in Eq. ([6]) is exponentially localized on 
low-energy levels, which acts as an effective high-energy 
cutoff [15(| . 

Although different modes are decoupled in the Ric- 
cati equations above, different z q (t) evolve with the same 
K(t). Therefore, noise-averaged many-mode quantities 
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must be computed by taking the noise-induced corre- 
lations into account. By expanding the above Riccati 
equation in 5K, we then obtain the following nonlinear 
Langevin equation: 



q 

K 



1) - q {K 2 z 2 + 1) 5a 



(7) 



where 8a(t) = —8K(t)/ Kg. Note that there is no dissi- 
pative term in the above equation. 

We now assume Gaussian noise with zero average and 
second moments characterized by strength W and corre- 
lation time t as in the Ornstcin-Uhlcnbeck process: 



8a(ti)5a{t2) 



2t 



D -\ti-t 2 \/-r 



(8) 



In the limit r — > , the right-hand side of Eq. (|8|) re- 
duced to W 2 8(ti — t 2 ) describing Gaussian white noise. 
We leave a discussion of the effects of colored noise (finite 
t) to the supplemental material [l5(, and only consider 
the r — > + limit here. As the white noise is thought of as 
a limit of a continuous process, we use the Stratonovich 
interpretation for Eq. ([7]) [u^ |. 

The stochastic differential equation described by 
Eqs. ((7j) and dU) is one of the key equations of this paper: 
it governs the stochastic evolution of z q (t), which in turn 
determines the many-body wave function [Eq. ([6])], and, 
consequently, all the observables of the system. We can 
write observables such as (H r q n (K)), where we have sup- 
pressed the K and 3 superscripts, in terms of z q by simple 
integration over the bosonic fields $ g . Similarly, equal- 
time correlation functions can be expressed in terms of 



15| . For example: 



i , 

q>0 



[q(x - x')] 
q$tz q 



(9) 



As we are interested in the low-energy limit of 
\Sz q (t)\ <C K q , where z q {t) = K^ 1 + Sz q (t), we can 
linearize Eq. ([7]) in 5z q and obtain a simple linear equa- 
tion 0: 

i 5z q = 2q (Szq — 5a) , (10) 

which admits the explicit solution 5z q (t) = 
2iq J* dt' e 2lq ^ ~ t " l 5a(t') (we will present an exact 
treatment of the nonlinear equation later). Our strategy 
for computing the noise average of generic functions 
of z q is then as follows. We expand these functions to 
leading order in 5z q , and insert the explicit solution 
above (in terms of 5a) into the resulting expressions. 
The averaging over noise can then be done by using the 
Wick's theorem, which relates 5a{ti)5a{t 2 ) ■ ■ ■ 5a(t n ) to 
two-point functions (|8]). (As the n-point function above 
vanishes for odd n, we may need to expand our function 
to one higher order in 5z q .) Finally, we perform the 



required integrations over the time arguments to find 
the noise average. 

Let us now compute the two types of energy moments. 
Here, we only consider the first and the second mo- 
ments (see the supplement for a discussion of higher mo- 
ments [15(). As stated above, all observables, including 
(ff™(-Ko)}> can be written in terms of z q . The expres- 
sions for (_ff™(-Ko)) become more and more involved as m 
increases, but if we expand these expressions to leading 
order in 5z q , we find the simple relationship (H^^Kq)) = 
q m 2 m ~ 3 K 2 \5z q \ 2 + 0{5zl). Using the strategy outlined 
after Eq. CEU]), we then obtain (H™) re 2 m - 1 q m+2 K$W 2 t 
and JWqY re q & K%W A [2t 2 + sin 2 (2gt)/4<j 2 ] [Hj]. Clearly, 
for a single mode, the energy fluctuations are significantly 
affected by the quantum source at least in the limit of 
small excess energy. 

Let us now consider the many-mode LL. The average 

H q ), where the 

factor of 2 accounts for the contributions of and H®. 



energy can be written as (e) = 2 X) 9 >o 

tioi 

Using J2 q >o -> ^JiT' wc then obtain 



(e) re Ltt 3 KlW 2 tj\. 



(11) 



We now relate the two types of second moments as fol- 
lows. We have (e 2 ) = (H 2 ) and (e) 2 = (H) 2 , with 
H = J2 q >o V^q These two moments can then 

be written, respectively, as sums over q\ and q 2 of terms 
of the form {H qi H q2 ) and (H qi ) (H q2 } , which only differ 
when qi = qi and both H qi and H q2 are in the same 5R 
or 3 sector (the wave function is a direct product). This 
leads to 



{e 2 ) = (e) 2 +2j2((H 2 q )~(H q ) 



(12) 



In the above expression, both (e) 2 and (e 2 ) scale as 
L 2 , while the sum over q scales as L. Therefore, if we 
take the thermodynamic limit before any other limit, the 
two types of moments will be, to leading order, identical. 
As we will see below, this is also the case for the two 



(0 



from 



types of cumulants obtained by subtracting 

the moments above. Such subtraction does not change 
the scaling with L 2 as an explicit calculation yields [15( : 



(e) 2 - (e) 



^KtwH 2 L 2 T{-Kt\ (13) 



where the function T[x) has the following asymptotic 
behavior: T{x) re 2 (l - \x 2 ) for and T{x) re ^ 

for i» 1. This indicates that the cumulant above crosses 
over from quadratic growth in t for short times to linear 
growth in t at longer times. Combining Eqs. (|13[) and 
(fTTj) leads to the important relationship ((T|). Recall that 
the regime of validity for these results can be extended 
by decreasing the strength of the noise W. 

We now turn to the effects of noise on the correlation 
functions. We compute the noise average of correlation 
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Correlation function Ground state Correction due to noise 
(<f>(x)<f>{x')) ln(s-s') (xW 2 t 3 (x-x')- 4 

{d x $(x)d x ,<i>(x')) (x-x')~ 2 ocW 2 t 3 (x-x')- 6 
(e^We-'"^')) \(x - x')-^/ 2 \(xW 2 t :i (x - x'Y^ I 2 

TABLE I: The leading scaling behavior of the noise-induced 
correction at time t to equal-time correlation functions. Stan- 
dard ground-state behavior [t = 0) is also shown. 



function (|9]), as well as those of current, d x Q(x), and ver- 
tex, e 4 ^*^), operators, which, for one realization of noise, 
can be derived from Eq. ((9]) by simple differentiation and 
exponentiation. We are interested in the asymptotic limit 
1 <C i - i' < I (in units of the lattice spacing). As the 
LL description is valid for finite t, we also assume 
Again, we use the method outlined after E q. (ITUl) . The 
calculations are tedious but straightforward [15[ , and the 
results for the scaling behavior of the corrections to these 
noise-averaged correlation functions are summarized in 
Table. |U Noting that a logarithm has scaling dimension 
zero, we find that these corrections exhibit the universal 
feature, Eq. (J2J. 

An alternative approach, which allows us to go beyond 
the limit of small 5z q , is through the Fokker-Planck (FP) 
equation for the wave-function probability distribution. 
Let us first briefly review the general formalism [Tt| : for 
a vector a of stochastic variables satisfying the Langevin 
equation d t a,i = hi(a) + g%{S)j(t), where hi and gt are 
arbitrary functions of a and 7(i)7(i') = 25 (t — t'), the 
probability distribution f(a,t) evolves according to the 
FP equation dtf = T>f, with the differential operator 
V = —sk^i ~ £-§^9 3 + £§7af-ftSj (summation over 
repeated indices is implied). The noise- averaged expecta- 
tion value of an arbitrary function G(a) of the stochastic 
variables a can then be computed at time t as an integral 
over J\i ddi weighted by the formal solution of the FP 
equation, f(a,t) = e vt f(a,0). 

In analogy with the Heisenberg picture of quantum 
dynamics, we can evolve the observable G(a) instead of 
the distribution function f(a) (using repeated integration 
by parts), and write 

W)\ t = I II M - f( S >°) ePtt G(a), (14) 

i 

where X* = (hi + f^-fifj) £- + 9i9j £: By expand- 
ing the exponential operator above as e ^ = ^2 n jjT>^ n , 
we can then compute the noise average of G(a) as a Tay- 
lor expansion in t. Since we are interested in finite time 
scales, such expansion is indeed very useful even when 
truncated at a finite order. 

For the average energy (e), it turns out, however, that 
we can resum the Taylor scries and obtain an exact so- 
lution. As the energy is just a sum of the single-mode 



energies, we need a FP equation for a single mode, i.e., 
a = (5Rz g ,32: 9 ). We can then explicitly construct the 
differential operator Tft corresponding to Eq. (J7J) [HI, 
and find the nth order term by applying it n times to 
the exact expression for (H q ), which can be easily writ- 
ten in terms of a [l5|. Note that the integration over 
a» in Eq. (|14p is trivial due to a delta-function initial 
distribution, /(a, 0) = 5(ai — Kq 1 )5((12)- By explicitly 
computing the terms in such expansion to order n = 10, 
we found that the series matches a simple exponential, 
(Hq) = | [exp (2q 2 W 2 Klt) - l] , term by term. Model- 
ing Eq. ([7]) by a sequence of random quenches [lH, [H| , 
we have also checked this expression numerically. Upon 
integration over g, we then obtain 

As expected, the above expression, up to first order in 
time, agrees with Eq. (fTT|) . 

Using the Wigner-function approach [3] (see also 
[ill ] ) , we can additionally show that for a thermal initial 
state (which is subsequently decoupled from the thermal 
environment during the evolution), the absorbed single- 
mode energy above is simply multiplied by a prefactor 
coih {q/ 2k bTq), where To is the initial temperature, and 
ks the Boltzmann constant [Hj]. The many- mode en- 
ergy can then be similarly computed by integration over 
q. Our FP approach is most useful for quantities such as 
the first and the second moments of energy and current- 
current correlation functions, where the noise-induced 
correlations do not need to be treated simultaneously for 
a large number of modes. The vertex-operator corre- 
lation functions, e.g., are not easily tractable with this 
method. Note that, to compute the second moment of 
the total energy, one needs to construct a FP equation 
for four stochastic variable, a = (^tz qi ,^sz qi ,^tzq 2 ,^sZq 2 ), 
and so on and so forth for higher moments. 

In summary, we studied a thermally isolated LL, ran- 
domly driven with a noisy Luttinger parameter, and un- 
dergoing coherent quantum evolution for each realiza- 
tion of noise. We computed noise-averaged correlation 
functions, and studied the energy fluctuations. We char- 
acterized these fluctuations by two types of energy mo- 
ments: one that mixes the classical and quantum sources 
of fluctuations and one that is only affected by the clas- 
sical source. We found that while for a single mode, 
the two types of moments lead to very different results, 
for the many-mode problem, the difference disappears 
in the thermodynamic limit. This indicates that many- 
body properties of such noise-driven coherent systems 
likely exhibit effective decoherence. Our approach to the 
dynamics of noisy LLs is based on a mapping of the quan- 
tum problem to one of nonequilibrium classical statistical 
mechanics, which provides powerful tools, such as the FP 
equation, for performing exact calculations. 
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I. EIGENSTATE EXPANSION OF THE SINGLE-MODE WAVE FUNCTION 

In this section, we compute the overlaps of the single-mode wave function with the eigenstates of the single-mode 
Hamiltonian H q . The wave function is given by 

i 

M4>, *) = ( 2g[5R / ?W] ) ' exp (-q z q (t) <f) , (1) 

1 where <b rcDrcscnts cither 



too 



where </> represents either <& q or <I>^ [see Eq. (6) of the main text]. Analogously, the Hamiltonian H q is shorthand for 



either or 



Q H In the absence of noise, the Hamiltonian is given by H = H q (K a ). For z q = 1/K , the wave function ([T]) is the ground 
state of Hamiltonian @. For arbitrary z q , however, it is a superposition ip q (4>) = Y^Lo c n(z q )ip n (4>) , where ip n (4>) is 
an eigenfunction of Ho with energy (n + ^) uq, which can be written explicitly in terms of the Hermite polynomials. 
Through direct integration, we can compute the amplitudes c n (z q ) = J dcf> ^p q (4>)*tjj n (<f)), which yield: 

I / \i2 r, I / M 2 (2n)! x\' 2 [{x q -lf + ^} n 

\c 2n+1 (z q )\ =0, =^TT^ [{Xq + 1)2+y2q] n + l/2> 

where 

x q = K [3? z g (i))] , y q = K Q [3 z q (t)] . 

We see from the above expression that the overlaps |c ra (z g )| are identically zero for odd n, and decay exponentially 
for even n. 

II. EXPRESSIONS FOR OBSERVABLES 

In this section, we express some observables in terms of z q (t). Measuring the energies with respect to the ground 
state of H q (Kq), and using the explicit form of the wave function ([1]), we can obtain (H q (Ko)) by direct integration 
over the bosonic fields: 



r+oo 



'l + A n 2 |zJ 2 ) -1 



2K^z q " 



(3) 



where u is set to unity. Similarly we can compute the second moment of energy: 



(H 2 q (K )) = (H q (K )) 2 + [1 - 2A 2 (5Rz 2 - 3z 2 ) + A 4 |z ? | 4 ] . (4) 

Note that both expressions above vanish for the ground state (z q = Kq 1 ). This simple approach can be used to 
compute higher moments of energy, but the resulting expressions become more cumbersome. Using MATHEMATICA®, 
we have calculated the moments (H q l (Ko)} for m < 20, and checked that to leading order in Sz q = z q — Kq , they 
can be written as: 

(O = (H™(K )) = q m 2 m ~ 3 K*\6z 9 \ 2 + 0(5z 3 q ). (5) 
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Equal-time correlation functions arc also simple to compute in term of z q (t). For example, Eq. (9) of 
the main text can be derived as follows. We can expand in Fourier modes, = J2 q ^9 = 

S g >o [cos(gx)3?$ 9 — sin(ga;)Sj$ g ], and make use of the relation (3fi$ gi 3fJ$ 92 ) = (3$ gi 3$ g2 ) = 8 qiq2 /Aq^.z q to 
write 

\ ^ / /^^ 1 V — > COS \q(x — x')] ,„s 

For the ground state, z 9 = , Eq. © reduces to its usual asymptotic logarithmic form in the thermodynamic limit 
XT*): 

($( a )*(0)> GS = g^dg^M = g(Ci( TO ) - Ci(ex)), (7) 

where e ~ ^ is a short-distance momentum cutoff, and the cosine integral Ci(C) = — ^j-^dt has the following 
asymptotic form for small £: Ci(C) ~ 7 + ln£, where 7 is the Eulcr-Mascheroni constant. In the asymptotic limit, 
1<C 1 C L, where the lattice spacing is set to unity, we can neglect Ci(7ra;), and recover the usual logarithmic 
correlation function. 



III. HIGHER MOMENTS OF ENERGY 



Focusing on a single mode for simplicity, we calculate the noise-averaged moments (e q ) m . To leading order in 8z qi 
we can write (e q ) m as [see Eq. above] 



\8z q \ 2 j . (8) 

After substituting the explicit solution, 8z q (t) — 2iq J Q dt' e 2tq( - t ~*-'i5a(i'), of the linear Langevin equation [see Eq. (10) 
of the main text] , we need to compute 



\8z q \ 2m = {2q) 2m f dt ie 2iqtl f dt[e- 2l < . . . f dt' m e~ 2i ^ x8a(t 1 )8a(t[)...8a(t m )8a(t' m ). 
Jo Jo Jo 



(9) 



For Gaussian white noise, 8a(t)8a(f) = W 2 8(t — t'), the above 2m-point function breaks up into a sum of products 
of to delta functions. The number of terms in this sum, i.e., the number of different contractions, is (2m — 1)!!. There 
are three different types of contractions: (i) 6a(ti)6a(tj), (ii) <5a(^)5a(^), and (hi) 6a(ti)8a(tj). Note that we have 
an equal number of contractions of type (i) and (ii), which are related to each other by complex conjugation. Upon 
integration, we obtain 



*dti f dtjeVto+VSaMSafo) x j* dt' m J* dt' n e- 2i ^ +t ^ 8a(t' m )8a(t' n ) = W 4 (^^) , 

t ft 

dU / dt'.e^'-^SaitiJSait'A = W 2 t. 



t 

to V 

J 

10 Jo 
We can then write 

Lm/2J 



(10) 



q 3m 



{KtwhT £ s m - p (^y\ (ii) 



p=0 



where [^J indicates the floor value of x and S m ,p represents the number of ways for having p contractions of type (i), 
p contractions of type (ii), and, consequently, m — 2p contractions of type (iii). The simple combinatorial argument 
below gives 



m 
2p 



(2p-l)!! 



2 



(to — 2p)\. 



Out of the to choices for Sa(ti), we select 2p, and then make p contractions among them in (2p— 1)!! ways. A similar 
argument applies to <5a(^), which results in the power of 2 for the term in the square brackets. We now have to — 2p 
unpaired 8a(ti), and to — 2p unpaired 8a.(t'j), which can be paired up in (to — 2p)\ ways. 
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IV. DERIVATION OF EQ. (13) OF THE MAIN TEXT 

Let us now explicitly compute (e) 2 . Since the energy is the sum of single-modes energies, (e) = 2^2 q (H q ) (the 
factor of 2 accounts for the contribution of H q and H q ), we can write (e) 2 = ^J2 qi ?2 {Hqi){Hq 2 ). After inserting 
Eq. (J5J for m = 1, we obtain 
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\Sz qi \ 2 \Sz q2 \ 2 , (12) 

91,92 

which, upon substituting the explicit solution of the linear Langevin equation [see Eq. (10) of the main text], yields 
W w 4A' 4 V q\ql [[[[ dhdhdt'^ e^^^-'^e^^-'i^a^Orfa^)^^)^^). (13) 

91,92 JJJJ ° 



As we did in in Sec. IIII1 we use the Wick's theorem to break the 4-point function Sa(ti)Sa(t' 1 )5a(t2)Sa(t' 2 ) into a 
sum of products of 2-points functions. However, in Sec. IIII1 we only had one mode q, and two types of contractions. 
Here, on the other hand, we have three different types of contractions. Upon inserting these delta functions and 
performing the integrals over time, we find 

W^j:j4{^^,^^ + e}. (i4) 



The term proportional to t 2 is exactly equal to (e) . To obtain Eq. (13) of the main text, we convert the sums to 
integrals, and define new variables qi = 7r£; and x = nt, which lead to Eq. (13) of the main text with 

m - 1 6 r «, t «, && I ±m + si °; [fe =m } . m 



JO 



For x <C 1, the asymptotic behavior is obtained by a simple Taylor expansion in x, which gives F(x) ~ 
2(l — |x 2 ). For x 3> 1, the main contribution comes from £i « £ 2 then we can approximate J-{x) ~ 

16 g J dQ J dq ^^j- 2 ^ (^ 2 ^) "'(gi)^^ wnere we have defined Q = £i + £ 2 an d <? = £i — £2 an d the factor 1/2 comes 
from the Jacobian of the transformation. Since Q 3> q the dominant contribution comes from the term proportional 
to Q 6 . The final change of variable k = qx leads to F(x) ~ ±£ J Q 2 dQ ^ QS s jf {k) = 



V. EFFECTS OF COLORED NOISE 

In this section, we discuss the effects of finite correlation time r in the Ornstein-Uhlenbeck process: 

<ya(ti)(5a(* / 1 ) = -_e-l 4l -*il/ T . (16) 



For the single-mode energy moments (e™) « q m 2 m 3 K 2 \Sz q \ 2 , we need to evaluate 

\8^={2q) 2 ( dh f dt^e 21 ^-^ 5a(h)Sa(t[). (17) 



(18) 



JO Jo 

Substituting Eq. (fTB"|). and performing the integral yields: 



{\Sz q \ 2 ) = Aq'WH 



1 r / e -T( 1 + 2 9T) _ 1 



1 + (2 9 t) 2 t \ (l + 2igr) 2 



Note that, for t>T, the above expression simplifies to |<5;z g | 2 = 4g 2 1+ ^ qT yi t, while for white noise, r — > + , we obtain 
\8z q \ 2 — Aq 2 W 2 t. We then conclude that, for t ^s> r, the single-mode energy can be obtained from the corresponding 
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expression for white noise by the simple rescaling W 2 — > 1+ ^ qT y^ ■ The average energy can then be calculated by 
integration over momentum: 



H = 2 IZ"S7~ [ « ?TT |^ = L ~ l^rf - log(l + Chrft] , 



(19) 



which for white noise, r — >■ + , reduces to (e) = Ln 3 'K^W^t / '4. 

Let us now consider the higher moments of the single-mode energy, which are given by (e q ) 
Eq. flSJ]. Using the method of Sec. IIII1 we arrive at the analogue of Eq. ([TO]) , which reads 



<ft { / dt 3 e 2iq{u+t ^8 a {U)5a{t. } ) x dt' m dt' n e- 2iq( - tm+u] ' 6a(t' m )Sa(t' n ) 



2qr + e'/ r (sin(2gt) - 2grcos(2gi)) 



^f\5z q \A [see 



W A e- 2t ' T 



2g(l + (2gr) 2 ) 
dti / dty^-^SaiQSaitr) = W 2 t 



(20) 



o Jo 



1 



■A(1+2«t) _ x 



l + (2gr) 2 i \ (l + 2^gT) 2 



Note that, for t ^> r, the two expressions above respectively simplify to 



W 2 



,l + (2gr) 2 
while for white noise, we obtain 



sin(2gi) — 2grcos(2gt) 
2^ 



W 2 



l + (2gr) 2 ' 



W 



sm(2qt)\ 



W 2 t 



(21) 



(22) 



by replacing r = in the above expressions. 

We observe that finite correlation time does not qualitatively change the white-noise behavior. We can now write 
the analogue of Eq. (fTTj) as 



W 2 



2 \ m L m / 2 J 



1 + (2gr)< 



p=0 



sin(2gi) — 2qTCOs(2qt) 
2qt 



2l> 



(23) 



where S m ,p is defined in Sec. IIIII 

Finally, we generalize Eq. (fT5|) . By repeating the derivation in Sec.|lVl we arrive at the analogue of Eq. (fl"4)) . which 
in the limit of t ^> r reduces to 



<e) 2 =4A-> 4 £ 



[i(l + ml) + T(m 2 - 1)] [t(l + m 2 ) + r(m 2 - 1)] 



(l + m 2 ) 2 (l + m 2 ) 2 



+ 



fa + ggMj- + miw 2 ) cos [fa + q 2 )t] - (l + ^(m 2 + m|)) sin [fa + g 2 )f] 
fa + 92 )(l + m 2 )(l + m 2 ) 

(gi ~ 92)t(1 + mim 2 ) cos [fa - g 2 )f] - (l + \{m\ + m\)) sin [fa - g 2 )t] 
fa -g 2 )(l + ™?)(l + m 2 ) 



+ 



(24) 



where we have introduced m, = 2g^r. The first term in the curly brackets is equal to 

see that in the limit r — > + we recover Eq. (HU). Replacing summation over g.; by 
variables: 



we obtain 



gi = 7r^, x = irt, y = 7tt, 



2 1 



(e) 2 = (e) + -L 2 K? ) WH 2 ^F{x 1 y), 



Moreover, it is easy to 
dqi, and making a change of 

(25) 
(26) 



5 



515 



2(1 



~ A ^I*** 4AA 



y = 0.0 

y = o.i 
y = 0.2 

y = 0.5 



\ 1 

12 3 4 

x 

FIG. 1: Eq. (|27[) for y = 7rr = 0, 0.1, 0.2, 0.5. Note that the equation is valid for x 2> y- The y = line corresponds to the 
white-noise case. The asymptotic behavior for i< 1 and y = is also shown (solid black line). 



where 

' ( (6 + 6Mi + 4ft6y 2 ) cos [(a + fejg] - (i + 2(aV + e 2 V)) sin [(a + j 2 )x] 
V ^i+6)(i + (26y) 2 )(i + (26y) 2 ) 

/ (6 - 6Mi - 466y 2 ) cos [(a - 6)*] - (i + 2(aV + £ 2 V )) sin [(a - 6K 

I x(a-a)(i + (2a2/) 2 )(i + (2ay) 2 ) 

Note that this expression is valid only for i>r, which implies The asymptotic behavior at long times, x 3> 1, 

can be obtained as for Eq. (fT5|): F(x,y) ~ (l — 4^y 2 + C(y 4 )). This implies that the finite correlation time, 
y > 0, does not change the form of the time dependence and simply lowers, by an overall coefficient, the white-noise 
(y = 0) asymptotic value. The short-time behavior can be seen in Fig. [T] for several values of y. In the physically 
relevant region, i>i/, the curves show a rapid decay followed by a slower one. This behavior is qualitatively similar 
to the white- noise case [Eq. (fl~5|) ]. which is also plotted in Fig. Q] for comparison. 





VI. NONLINEAR FOKKER-PLANCK EQUATION 



In this section, we explicitly derive the differential operator Tfi used in the derivation of Eq. (15) of the main text 
[see also Eq. (14) of the main text]. Let us define 

We can then write the nonlinear Langevin equation [Eq. (7) of the main text] as 

&i = 2K q t MiA - 2Kl qi MiJi 8a, 

A = K oqi {J 2 - ®\ + AV 2 ) - Kfa {J 2 - *? AV 2 ) 5a. 

As the total energy is a sum of single-mode energies, we only need a single-mode Fokkcr-Planck equation, and we can 
drop the subscript i. Rescaling 5a —> ^7(*)i with (j{t)j(t')) = 26(t — t'), we can read off /j and Qi for i = 1,2 (in 
the notation of the main text): 

f 1 =2K qa 1 a 2 , gi = -V2W K 2 q ai a 2 , f 2 = K q (a 2 - a 2 + K^ 2 ) , g 2 = -^=K 2 q (a 2 - aj - K^ 2 ) , (29) 

where a= (ai,a 2 ) = [ffl,^). 
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In terms of the above /, and Qi , we can write 



=D[ 1 >a ai +Dl l >d«+D['{% 1 



where 



D 



- /l + .9i ~ h.92 



9gi 

! "5 — 

oa 2 



D 



(i) 



r>(2) S 2 

LJ 2.2 u a 2 



dg2 



f2 + gi + 92 



dg2 

! "5 — 

oa 2 



D 



(2) 



mi- 



(30) 



(31) 



We now consider the single-mode energy [Eq. ([3])]. The general expression given in the main text [see Eq. (12)] gives 
the Taylor scries of the noise-averaged quantity: (e q ) = X) n =o flT c ™' wnere the coefficient c„ is obtained by applying 
the operator T>^ n to Eq. ([3]), and setting 5Rz — Kq 1 and 3z = in the resulting expression. We have explicitly 
calculated such terms for n < 10, and found that cq = 0, and, for n > 1, c n = ^ (2q 2 W 2 K^) n . We are thus lead 
to conjecture that the behavior holds for all n, and the series can be resummed to the exponential expression in the 
main text. We have also checked this conjecture by numerical simulations. Note that if we are interested in 2-modcs 
quantities, we need to repeat the procedure above with a = &2, ^2), identify fi,gi for i = 1,2, 3, 4, and build 

the 4-dimensional vector and the 4x4 matrix D^ 2 \ In general, if we need the correlations of m momenta, 
is 2m-dimensional and is a 2m x 2m matrix. 



VII. WIGNER-FUNCTION APPROACH 



An alternative approach to the dynamics of the problem is through the Wigner-function representation. This 
approach allows us to treat, on the same footing, initial conditions given by a finite-temperature thermal density 
matrix as well as the zero-temperature ground state considered so far. Note that the system is still thermally isolated 
during the evolution, i.e., it is decoupled from a finitc-tcmpcrature heat bath at time t = 0. 

Let us briefly review the formalism. In terms of phase-space variable x and p (which can be vectors for multi- 
dimensional problems), the Weyl symbol £l w (x,p) for a quantum operator O is defined as 



Q w {x,p) 









i 


I ds ( x 

/ \ 2 


n 


x+ CX P 





(32) 



The Wigner function W(x,p) is, by definition, the Weyl symbol for the density matrix p, which satisfies the following 
equation of motion: 



d t W = ~H w sm(^A] 



(33) 



where H w is the Weyl symbol for the Hamiltonian, and A = dp. The expectation value of an operator O 

can be written as 



(6> = 



dxdp 
2tt 



W(x,p) O w {x,p) 



(34) 



in terms of its Weyl symbol O w {x,p) and the Wigner function W(x,p). 

Let us now consider the single-mode Hamiltonian H q (K(t)) = (j^p- p 2 + j^pj q 2 x 2 ^j [sec Eq. ©], where we have 

set u = 1 and suppressed the 3? or 3 superscripts. Here, x and p respectively represent the real (or imaginary) part of 
operators <fr q and H q , and we have used the hat notation to tell apart quantum operators from phase-space variables. 
For a system evolving from an initial thermal density matrix, 

(35) 



(36) 



Pa = cxp \-H q {Ko)/k B T j fix [cxp (-H q (K )/k B T 
with the Hamiltonian above (for arbitrary K(t)), the Wigner function retains the following form: 

W(x,p) = Afexp 

where A, B, C and M arc potentially time-dependent functions with the following initial conditions 



A 2 B 2 
-—x - — p +Cxp 



A(t = 0) = 



KqR' 



B(t = 0) = ^ C(t = 0) = 0, jV(t = 0) = |, H=coth 



(37) 
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Note that the normalization condition JJ^^ ^^-W(x,p) — 1 sets 



TV = y/AB - C 2 (38) 

at all times. For the quadratic hamiltonian H q (K(t)), the Weyl symbol is simply obtained by replacing the operators 
x,p with the phase-space variable x,p. We can then write Eq. (|33p as 



Af A^ 2 B 2 nrt --\( n I 2 -in. D _1 ^(*) 



JU(x,p) I ^ - -x z - -v* + Cxpj = W(x,p) 1^2^ x (Cx -Bp)--^p (Cp - Ax) 

which leads to the following equations for the parameters of the Wigner function: 

tf = 0, A = -4^-C, B = K(t)C, C = 2{^A-^ )B ). (39) 

By using Eq. ()38|) , and noting that Af does not change from its initial value, we can eliminate A from the equations 
above: 

We now expand Eq. (|40[) in SK(t) <C Kq, and obtain two coupled Langevin equations: 

B = CK - CK$5a(t), 



C = 



K0(4 + C 2 R 2 ) 2q 2 B 
2R 2 B 



2Bq2+ K 2 { , + C 2 R 2 ) 



2BR 2 



5a(t), 



(41) 



where 5a(t) = —5K(t)/K 2 . Using Eqs. (|34j) and ([37| , we can also express the absorbed energy in terms B and C as 

- , _ K (A + C 2 R 2 ) q 2 R 2 B qR 
(H q ) (H q ) t = — + - — . (42) 

As expected, for R = 1 (T = 0), Eqs. (|4"Tj) and (l4"2"j) above are respectively equivalent to Eqs. (J2HJ) and ^ through a 
change of variables: B = l/q$tz q and C = —2^sz q /diz q . Applying the Fokker-Planck approach of Sec. IVII to Eqs. (|4Tj) 
and (j42|) . we obtain (e q ) = ^ (cxp [2fcQpU 2 q 2 i] — l), and find that the noise-averaged absorbed energy depends on 
the initial temperature only through the prcfactor R = coth ^^fefrb) ' 

VIII. NOISE-AVERAGED CORRELATION FUNCTIONS 

In this section, working in the limit of small 5z q , we consider the effect of noise on the following noise-averaged 
correlation functions: 



d(x-x') = (®(x)$(x')), C 2 (x-x') = (d x $(x) d x ,<f>(x')), C 3 {x - x') = {e iv ^ x ) e~ iv ^ x ')) . (43) 
Notice that even though, in the absence of noise (or for each realization of noise), the three correlation functions above 



are related by simple differentiation or exponentiation, i.e., C 2 (x — x') = d x d x i C\{x — x') and C^(x — x') = e~ Cl ( x ~ x ' , 
the noise averaging must be done at the last step, and such simple relations do not necessarily hold for noise-averaged 
quantities. 

Let us start with C\(x) = jr J2 q>0 c °^ x ^ [see Eq. ©]. Bearing in mind that the noise average of an odd power of 
6z q vanishes, to leading order in 8z q , we can write 



C i (*) ~ j; U K ^1L ( 1 + K 2 (mz q f ) . (44) 



1 ^--^ T ^ cos(qx) 

q>0 

We can now write (mz q ) 2 = 4q 2 W 2 f * sin 2 [2q(f - £)] dt' = ^-q [Aqt - sm{Aqt)\, which yields 



d(x) « (K J"dq^& + ^1 £ dqcos(qx) [4gt - sin(4gi)]) , 
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where e ~ ^ is a short-distance momentum cutoff. The first term in the right-hand side of the above equation is the 
correlation function c[°\x) in the ground state [see Eq. (|7|)]. We can then perform the integral in the second term 



and write 



(o) K§W 2 f 8t - (At + x) cos[7r(4i - x)} - (At - x) cos[7r(4i + x)] 4t[cos(7ra;) + 7rx sm(7rx) - 1] 



4tt V 2a; 2 - 32i 2 



(46) 



The right-hand side of the above equation simplifies by neglecting fast oscillatory terms for i>1, and using i»t 
(small Sz q requires finite t): 

Notice that we could obtain the same result, without an exact computation, by noting that, in the asymptotic limit of 
large x, the main contribution comes from small momenta. We can then expand Aqt — sin(4q<) ~ (Aqt) 3 /3 in Eq. (|4"5")) 
for small q, and compute a much simpler integral, which directly leads to Eq. (|47[) . 
Similarly, we can write 

C^x) « i I<0 1 cos (9' T ) i 1 + K o P^F) - (48) 

(J>0 



which leads to 

i / r v K 3 W 2 f n 

dq q 2 cos(qx) [Aqt - sm(Aqt)} ) . (49) 



C z( x ) ~ ^~ (^ K o J dqqcos(qx) 



2 

The first integral above gives Q (x) = Ip 1 [—1 + cos(7ra;) + ttx sin(7rx)] /x 2 , which after neglecting the fast oscillatory 
terms reduced to the usual — scaling of the current-current correlation functions. By performing the second 

integral, and similarly neglecting the fast oscillatory terms, we can write the noise-induced correction as 



4tt V x4 (x + At) 3 (x-At) 3 



7T X" 



The calculation of C3 (x) is also possible by exactly the same method. The first step is to expand the correlation 
function to second order in 8z q . We begin by writing 



which gives 

C^x)^cf\x)Y[ 



q>0 



(1 



1 + ^rifg— ^ (-K<5z g + A (9Wz 9 ) 2 ) + -^K$ $-L (mz q ) 2 



. (52) 

~C$\x)\l + ^KlY,^W^+^Kt £ COS(glX)C0S(g2X) (^)(^) ■ 

I 9>0 9 ?i,«a>0 9192 J 

Notice that to go from the first to the second line, we have neglected a term of the form S 9 >o (which scales as 
-^). We can also write 



(mz qi )(mz q2 ) =A qi q 2 W 2 f sin[2gi(t'-t)]sin[2g 2 (t'-t)]di' = qmW 2 

Jo 



sin[2(q 1 - q 2 )t] sin[2(gi + q 2 )t] 



qi - 92 9i + 92 



(53) 

For large x, we can expand the above expression in small momenta, and obtain (#t5z qi )($l5z q2 ) ~ 3 ^-W 2 q 2 q 2 t 3 , which 
upon performing the integrals gives 

Cjx)-C™ (x)~t 3 /x^ 2 / 2 7 

with cf\x) = etrS^oft^ ^ . 



